Kinks and Mid-Infrared Optical Conductivity from Strong Electron Correlation 
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We compute the one-particle spectral function and the optical conductivity for the 2-d Hubbard 
model on a square lattice. The computational method is cellular dynamical mean- field theory 
(CDMFT) in which a 4-site Hubbard plaquette is embedded in a self- consistent bath. We obtain a 
'kink' feature in the dispersion of the spectral function and a mid-infrared (mid-IR) absorption peak 
in the optical conductivity, consistent with experimental data. Of the 256 plaquette states, only a 
single state which has d x 2_ y 2 symmetry contributes to the mid-IR, thereby suggesting a direct link 
with the pseudogap. Local correlations between doubly and singly occupied sites which lower the 
kinetic energy of a hole are the efficient cause of this effect. 



A distinctive feature seen in angle-resolved photoe- 
mission spectroscopy (ARPES) spectra on the high- 
temperature copper-oxide superconductors (cuprates) is 
the 'kink' 1 or abrupt change in the slope of the electron 
dispersion uo vs k. The abrupt change in the slope in 
the electron dispersion or the quasiparticle velocity ob- 
tained from ARPES, has been observed in various fam- 
ilies of hole-doped cuprates at an energy of 50 meV - 
80 meV. The size of the kink decreases with hole dop- 
ing x, and the kink seems to disappear around x = 0.3. 
Since ARPES offers a direct measure of the spectral func- 
tion, the imaginary part of the self-energy can also be 
deduced. Experimentally, the imaginary part of the self 
energy IrriE(k,w) shows a suppression, indicating a di- 
minished scattering rate 2 , at the kink energy. The pres- 
ence of the kink is reminiscent of similar observations in 
the spectra of conventional BCS superconductors, where 
phonons cause an abrupt change in the electron disper- 
sion as well. Consequently, explanations based on col- 
lective phonon excitations^ 3 -^^^ have been invoked 
to explain the ubiquitous kink- feature. Rival explana- 
tions include purely electronic scenarios based on spin- 
fluctuations 9 -^. Recent ARPES experiment s 11 ! 12 find 
that kink phenomena occur at high energies as well. At 
the high-energy kink (referred to colloquially as the wa- 
terfall), the electron dispersion bifurcates. 

Equally striking is the optical conductivity in the mid- 
infrared. The optical conductivity of a material can be 
extracted from its measured reflectivity R(uo) spectrum 
using the Kramers- Kronig (KK) relation s 13 ! 14 ! 19 . Two 
universal features seen in the optical conductivity of the 
cuprates are: l)a Drude peak at low energies followed 
by 2) an absorption feature in the mid-IR region. The 
Drude peak corresponds to the motion of free carriers in 
the system for which the Drude model predicts a con- 
ductivity peak at zero frequency. The existence of an 
absorption feature in the mid-IR region is unexpected, 
since doped Mott insulators are expected to have spectral 
weight either at the high energy sector across the Mott 
gap(UV) or at low energy, close to the Fermi level (far 
IR), 15 . The spectral weight of the mid-IR band in the 
optical conductivity is an increasing 16 ,17 ! 18 ! 20 ! 21 ! 22 func- 
tion of doping whereas the band at high energy (UV) 
decreases. This tradeoff suggests that the origin of the 



mid-IR band arises fundamentally from the strong cor- 
relations that mediate the Mott state. It is not sur- 
prising then that no low T c materials exhibit a mid-IR 
band (MIB) in their optical conductivity. Amidst the 
myriad explanations propose d 23 ! 24 ! 25 ! 26 ! 27 ! 28 ! 29 , none has 
emerged as the definitive answer to this problem. 

In a series of recent paper s 33 ! 34 ! 35 we formulated the 
exact low-energy theory of a doped Mott insulator and 
showed that a new excitation emerges at low energies. 
The excitation is a charge 2e boson which is not made out 
of elemental excitations. This excitation was shown to of- 
fer a unifying mechanism for the kink features (both high 
and low energy as well as the bifurcation of the electron 
dispersion) and the mid-IR band. What we present here 
are a series of numerical studies on the Hubbard model 
which offers an independent test of the predictions of the 
low-energy theory. To this end, we employ the state-of- 
the-art numerical method on the Hubbard model, cellu- 
lar dynamical mean- field theory (CDMFT) 30 , to investi- 
gate whether the low-energy kink and mid-IR bands have 
purely electronic origins. We find that both the kink and 
the mid-IR are linked to short-range electronic correla- 
tions. In particular, the mid-IR arises from a mixing 
of the upper and lower Hubbard bands through an elec- 
tronic state that has d x 2_ y 2 symmetry, suggesting that 
the pseudogap is also due to the same mechanism. This 
mechanism is identical to that mediated by the charge 2e 
boson in the exact low-energy theory. 

To this end, we start with the 2D Hubbard model, 
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. Here z, j label lattice sites, is equal to one if z, j 
are nearest neighbours, Ci a annihilates an electron with 
spin a on lattice site z, t is the nearest-neighbour hop- 
ping matrix element and U the energy cost when two 
electrons doubly occupy the same site. The cuprates live 
in the strongly coupled regime in which the interactions 
dominate as t ~ 0.5eV and U = 4eV. Various numerical 
techniques have been developed to study strongly inter- 
acting systems on a lattice. These include exact diago- 
nalization (Lanczos technique) , Quantum Monte Carlo 
and cluster methods. In this study, we use Cluster Dy- 
namical Mean Field Theory (CDMFT) 31 to compute the 
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FIG. 1: Intensity plot of the spectral function computed along 
the nodal r = (0, 0) -> Y = (tt,tt) direction for U = 8t at 
three different hole dopings. The arrows in the first two panels 
indicate the position of the kink. 



one-particle Green function. In this method, we consider 
exactly the dynamics on a 4-site (plaquette) cluster and 
treat the rest of the lattice as a bath whose degrees of 
freedom are integrated over. The coupling of the cluster 
to the bath is thus treated in a mean-field fashion and 
the cluster quantities are evaluated self consistently us- 
ing the cumulant lattice reconstruction scheme 3 -^. While 
no variational principle exists for this method, several 
tests in limiting cases in which the answers are known 
exactly have revealed that the local correlations that the 
CDMFT method builds in provide an adequate starting 
point for quantifying the physics of strong correlations. 

Using the CDMFT method, we obtained the one- 
particle spectral function A(k, uS) for various values of 
hole doping x and on-site Coulomb repulsion U (ex- 
pressed in units of the nearest neighbor hopping integral 
t) . Hole dopings in the range of 0.03 (lightly doped Mott 
insulator) up to around 0.30 (heavily overdoped super- 
conductor) have been considered together with a maxi- 
mum U value of 12t. Typical spectral functions, A(k, a;), 
obtained from CDMFT are displayed in three panels in 
Fig. [TJ We have plotted here only the spectral function 
intensity for different hole dopings in which the bright- 
ness at each point (k,o;) is indicative of the magnitude 
of A(k, uj) at that point. The bright band in each panel 
is the region of maximum A(k, u) and defines the uj vs 
k dispersion curve. The most interesting feature in these 
plots is the presence of a distinct kink in the dispersion 
curve at low doping values and an absence of the kink 
at the highest doping. The overall doping dependence is 
in rough agreement with the experimental trends. Also 
shown in Fig. [2] is the U-dependence of the kink energy 
for three U = 4t , U = 8t and U = 12t. As seen in 
Fig. [21 the kink energy scales inversely with U . Fitting 
the kink energies to 1/U gives a very good linear plot 
with as slope approximately 4. In other words the kink 
energy scale is given by At 2 /U . 

To explore the physics of the kink in further detail, 



FIG. 2: Left panel shows dispersion curves at x = 0.05 plotted 
for 3 different U values with the black arrows showing the 
position of the kink in each curve. Right panel shows kink 
energy plotted vs 1/U and its linear fit. 



we looked for specific short-range electron correlations 
within the plaquette that give rise to the kink. Such 
details are readily available from the 176 local resol- 
vents that comprise the impurity solver method we em- 
ployed, namely the non-crossing approximation (NCA). 
To achieve this, NCA has been selectively employed on 
various subspaces of the full Hilbert space of the plaque- 
tte in order to isolate the effect of each subspace on the 
spectral function. This requires an extensive search on 
the 4 4 = 256 dimensional Hilbert space and isolating the 
states relevant to kink formation. It turns out that out 
of 256 plaquette states, only 16 are involved in giving rise 
to the kink. All of these states ( labelled 'Super 16' in 
Fig. [3]) have 4 electrons with a total spin S z = ±1 i.e. 
there are three same-spin electrons and an opposite spin 
electron on a plaquette. All of the states in the 'Super 
16' subspace have some doubly occupied character. By 
doubly occupied character, we simply mean that there is 
a non-zero overlap with the doubly occupied sector not 
that the wavefunction is predominantly doubly occupied. 
We also verified that the kink vanishes if the doubly oc- 
cupied sector is eliminated. Consequently, the kink we 
have found here does not have a simple interpretation 
in terms of spin fluctuations as such a scenario would 
not require the explicitly doubly occupied sector. A kink 
arising from an explicit mixing between singly and doubly 
occupied sectors with S = 1 is consistent, however, with 
the physics mediated by the charge 2e boso n 33 i 34 i 35 that 
appears in the exact low-energy theory of a doped Mott 
insulator modeled by the Hubbard model. Excitation of 
the boson, which mediates mixing with all the doubly oc- 
cupied sectors, occurs on the energy scale t 2 /U. Once the 
boson is excited, the electron dispersion should change. 
Hence, we conclude that our numerical results are consis- 
tent with the physics that the charge 2e-boso n 33 i 35 me- 
diates at low energy in the Hubbard model. 

We also computed the optical conductivity. To obtain 
a direct link between the conductivity and the spectral 
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FIG. 3: Spectral function for x = 0.14 obtained in two re- 
stricted Hilbert spaces on the plaquette (center and right pan- 
els) compared with the complete spectrum (left panel). The 
kink goes away only when the 'Super 16' states are removed 
from the Hilbert space and reappears when they are included 
even when other states with S z =± 1 are absent. 



ing the artificial U = oo constraint results in the mixing 
of the high and low energy scale by virtue of a hole neigh- 
bouring a doubly occupied site. This resonance persists 
even at x = (though with perhaps vanishing weight) 
and hence we predict that the restricted f-sum rule at 
finite doping should extrapolate to a non-zero value at 
x = as long as U / oo as depicted in the inset of 
Fig. ([5]). As this is the only process through 0(t 2 /U) 
that is left out in the U = oo limit, we conclude that 
the resonance shown in Fig. ([6]) is the origin of the mid- 
IR. This would indicate that the mid-IR is strongly mo- 
mentum dependent (and hence not adequately described 
by single-site analyses 41 ), having d x 2_ y 2 symmetry, our 
principal conclusion. That the process in Fig. (j6]) causes 
the mid-IR was also the result of our analysis^ of the 
exact low-energy theory of the Hubbard model. 



function, we worked in the non-crossing approximation 

k 

xik^'+o;^ ^ J - (2) 

to the Kubo formula for the conductivit y 36 ! 37 ! 38 where 
f(uo) is the Fermi distribution function and A(lj, k) is the 
spectral function. Here v is the unit cell volume and ek is 
the dispersion for the non - interacting system. The op- 
tical conductivity has been computed for various dopings 
and U values. Fig. [H displays the optical conductivity 
for three different hole dopings. The optical conductiv- 
ity plots show a peak-like resonance feature (Fig. [4j) at 
0.5 eV (4000 cm -1 ) for U = St which falls right inside 
the mid- IR region where the experimental— i 17 i 18 i 2Q i 21 i 22 
data shows an absorption peak. The magnitude of the 
optical conductivity is also in perfect agreement with 
experimen t 17 ! 18 . The physical origin of this peak can be 
determined by focusing (see Fig. J5])) on the resolvents 
for states in the plaquette that contribute significantly 
to the optical conductivity. Surprisingly, of the 256 pla- 
quette states, only a single state in the N=4 and S z = 
sector contributes (see Fig. (|S])) to the mid-IR peak. This 
state has three key characteristics: 1) it contains a mix- 
ture of singly (87%) and doubly occupied (13%) sites, 2) 
its energy is — 1.3t, essentially the energy of mid-IR peak, 
and 3) the spatial symmetry of the eigenstate is d x 2_ y 2. 
Any physical process which meets these constraints must 
also couple to the charge not simply to the spin sector 
as in the case of antiferromagnetic spin fluctuations. A 
further hint as to the origin of the MIB is the calculation 
of Haule and Kotliar 39 indicating that the MIB is ab- 
sent in the traditional implementation of the t — J model 
unless superconductivity is present. By traditional im- 
plementation, we mean the U = oo limit in which J is 
(inconsistently) treated as being finite but double occu- 
pancy is excluded in the Hubbard basis rather than being 
excluded only from the transformed basis 4 ^. The typical 
process, shown in Fig. (|6|), which is eliminated by impos- 
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FIG. 4: (Color online). Optical conductivity plots computed 
using CDMFT for three different hole dopings for U = 8t, 
T = O.lt. Inset: Integrated optical conductivity K(Q) = 
J Q cr(uj)duj) up to energy of (Q = 60t = 30eV) (blue dots) 
and electron kinetic energy (red dots) plotted for three differ- 
ent hole dopings for U = 8t, T = O.lt. 

An explicit assumption in the NCA procedure (Eq. [2]) 
is that vertex corrections are not important. That is, the 
current-current correlator is approximated by the simple 
bubble diagram and higher order terms (vertex correc- 
tions) are neglected. This approximation is only correct 
in the limit of infinite dimensions^ 2 -. To determine how 
valid this approximation is for a d=2 system, we com- 
puted the kinetic energy explicitly versus the value pre- 
dicted from the sum rule 

(oo) = ~Wv Km { ) 

where K(Q) = a(co)duj is the integrated optical spec- 
tral weight, e is the electron charge, a is the lattice con- 
stant and V is the volume of the unit cell. This sum 
rule is expected to hold as long as the hopping terms are 
mediated by nearest-neighbour hopping processes. The 
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FIG. 5: Optical conductivity plot for U = St, n = 0.94, 
T = O.lt with distinct mid-IR peak contrasted with the 
contribution from a restricted Hilbert space with certain 
S=0 states removed. Inset: Integrated optical conductivity 
K(Q) — J a(uj)duj, up to (Q = 3t = l.beV) (restricted f- 
sum rule) for U — 12t and U = St. 



FIG. 6: Hopping process that mixes the upper and lower 
Hubbard bands that leads to the mid-IR band. This process 
requires that the hole and the doubly occupied site are neigh- 
bours. It is the motion of this bound state that we attribute 
to the mid-IR band. 



low hole dopings (x = 0.05) near half-filling. However, 
further away from half-filling where the low-energy sec- 
tor contains increased spectral weight transferred from 
high energy, the integrated optical conductivity fails to 
account for the kinetic energy. This behaviour is not un- 
expected as the vertex corrections ignored in Eq. ([2|) are 
strongly doping dependent and hence Eq. ([2]) is more 
accurate close to half-filling. As our focus is on the the 
mid-infrared feature which is most well-defined close to 
half-filling, our treatment here should be adequate. 

The key feature that distinguishes this work is the find- 
ing that double occupancy mediates both the low-energy 
kink and the mid-IR band. Our result that the mid-IR 
is caused by a single state (of the 256 plaquette states) 
that has d x i_ y i symmetry is striking and indicates that 
the mid-IR is likely to be caused by the same physics 
that underlies the pseudogap as optics experiments seem 
to indicate 14 . Further, the mechanism identified here, 
which is identical to that mediated by the the charge 
2e boson in the exact low-energy theory of the Hubbard 
mode l 33 i 34 i 35 (Fig. (|6j)) could explain why an extrapo- 
lated non-zero intercept of the restricted f-sum rule might 
be a generic feature of doped Mott insulators, in contrast 
to recent claims 41 . 
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